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ABSTRACT 

We present the leading order non-linear density and velocity power spectra in the com¬ 
plete form; previous studies have omitted the vector- and tensor-type perturbations simultane¬ 
ously excited by the scalar-type perturbation in nonlinear order. These additional contributions 
are comparable to the scalar-type purely relativistic perturbations, and thus negligible in the 
current paradigm of concordance cosmology: i.e., concerning density and velocity perturba¬ 
tions of the pressureless matter in perturbation regime well inside of matter-dominated epoch, 
we show that pure Einstein’s gravity contributions appearing from the third order are entirely 
negligible (five orders of magnitude smaller than the Newtonian contributions) in all scales. 
We thus prove that Newtonian perturbation theory is quite reliable in calculating the ampli¬ 
tude of matter fluctuations even in the precision era of cosmology. Therefore, besides the ones 
imprinted as the initial condition generated in the earlier phase, the other relativistic effect 
relevant for interpreting observational data must be the projection effect that occurs when 
mapping galaxies on to the observed coordinate. 

Key words: gravitation - hydrodynamics - relativity - cosmology: theory - large-scale struc¬ 
ture of Universe. 


1 INTRODUCTION 

Perturbation theory is a common practice in wide variety of physics problems. While nearly all fundamental equations in theoretical physics 
are non-linear, in most cases only the linear perturbation can be mathematically handled with sufficient confidence and precision. Thus, it 
is truly a surprise to discover lately that many of the cosmologically relevant observations can be accommodated in a simplest geometric 
background world model with additional linear perturbations. These are temperature anisotropies and polarizations in cosmic microwave 
background radiation and the distribution of galaxies in the large scale (large-scale structure) in the observational side, and the Friedmann 
world model with small perturbations in the theoretical side ( |Friedmann|| 1922 [ |Lifshitz|[l946| >, now encapsulated in the concordance or 
precision cosmology. Although their large-scale effects are supposed to be small in the concordance model, in this precision era of cosmology, 
the non-linear perturbations might be important so that we need to include higher order corrections. In fact, the fully nonlinear situation 
appears as we pay attention to the smaller scales in the large-scale structure. 

In cosmology, the linear perturbation is often decomposed into three types of perturbation. These are scalar-, vector- and tensor-types 
of perturbation jLifshitz|1946[[Wrk|1973[|Bardeen|1980[[ifodama & Sasaki|1984^ where the last one is absent in the Newtonian approach. 
Such a three-type decomposition, however, is mathematically not unique to the non-linear order jHwang & Noh|2013^ . Furthermore, as all 
three types of perturbations are convoluted in the equation level, the physical meaning of each type naturally becomes obscured as well. Only 
to the linear order and in the spatially homogeneous and isotropic background, the three types of perturbations decouple in the equation level 
( |Lifshitz|1946j . Thus, only to the linear-order perturbation in the Friedmann background world model the scalar-, vector- and tensor-type 
perturbations can be properly interpreted as corresponding to the density, the rotation, and the gravitational wave perturbations, respectively. 
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Even when we consider pure scalar-type perturbations to linear order, the vector- and tensor-type perturbations are naturally excited at 
non-linear order sourced by the linear scalar-type perturbation. 

The leading non-linear power spectra demands perturbation expanded to the third order ( |Vishniac]|1983^ . In our previous work, we 
presented the leading nonlinear power spectra of the density and velocity in the zero-pressure medium ( [Jeong et al.|2011^ . We considered 
Einstein’s gravity in the comoving gauge. In that work, we considered purely scalar-type perturbation, thus omitting the naturally accompa¬ 
nying effects of the vector- and tensor-type perturbations excited by the scalar-type one to the nonlinear order. 

Here, we include these effects which start appearing from the third-order in perturbation. We shall show that the second-order vector- 
and tensor-perturbations generated from the linear scalar-perturbations contribute to the third-order solutions of density and velocity pertur¬ 
bations. We derive the general relativistic energy and momentum conservation equations including these effects, see equations [16[ l and l |17[ l. 
In these equations, Y and Zij terms are the contributions from the vector- and tensor-perturbations, respectively. We present the complete 
leading non-linear power spectra of the density and velocity in Figs[2and|^ respectively. 


2 THIRD-ORDER PERTURBATIONS IN THE COMOVING GAUGE 


We consider zero-pressure fluid in a flat background, thus set p = 0 and K = 0 (variables are defined in Appendix]^. Our aim is to derive 
third-order perturbation equations for the density and velocity perturbations in the comoving gauge without linear-order vector and tensor 
perturbations. Thus, in what follows, we consider only scalar-type perturbations in linear order and vector- and tensor-type perturbations 
shall be treated as higher order quantities. We decompose the perturbation variables to the scalar and vector types as 


Vi = -Vi + V. 


■{") 


Xi = cX,i + X, 


(u) 


( 1 ) 


with = 0 = for our metric and energy-momentum convention, see Appendix [ a| By considering only scalar-type perturbations 

in the linear order, we set = 0 = to the linear order. From the second order, the vector-type perturbations are generated from the 
scalar perturbations. 

In our fully non-linear formulation jHwang & Noh|2013) , we have ignored the tensor-type perturbations. Just like the vector-type per¬ 
turbations are generated from the scalar perturbations in nonlinear order, the tensor-type perturbations must be generated as well. The tensor 
perturbation generated from the scalar ones is at least second order in perturbation. Thus, in our calculation, aiming for third-order solutions 
for density and velocity perturbations, we only need to keep the tensor-type perturbations appearing with the scalar-type perturbation: we 
will recover this contribution in equation below. Note that this is not the case for the vector perturbation as we shall show in what follows. 

By taking the comoving gauge, we set 


w = 0, 


( 2 ) 


to all orders in perturbation. 

We have identified S and k as the density and velocity perturbation variables, see equation l |l0| (; we have Jl = qc^, g = g + Sg and 
5g = gS. Equations for <5 and k follow from equations l |A9[ and ( |A7[ l, respectively. We note that the vector-type perturbations appear in the 
Vi and Xii in equations ( |A9^ and (|A7j, these terms appear at least in quadratic combination. As the vector perturbations generated from 
the scalar ones are at least second order, we need Vi and the vector part of Xi only to the second order. These can be read from equations 
| |A1Q| ( and jA6^ , respectively. 

For Vi, from equation jAlOj ), we have 

„2 




+ —Vi'^^ + —SAfi — non — linear terms all involving v 
a a ’ 


(v) 


( 3 ) 


to the fully non-linear order. The equation to the linear order reduces to 5Af = 0. To the second order, terms in the right-hand side vanish, 
thus we have &J\f = 0 and oc a~^. As^”^ is purely decaying proportional to the reciprocal of the scale factor in an expanding medium, 
we may set nl — 0 to the second order. This procedure can be continued perturbatively to all higher order perturbations, and we have 

5N = 0, = 0, thus Vi = 0, 


to the fully non-linear order. The vector part of Xi will be handled below. 
Equations l|A9|l and l|A7|( to the third order, respectively. 


5 — K H—^ (1 — ‘2‘fi) 5 iX — Sk = 0, 

k + 2 Hk — AirGgS + -^ (1 — 2ifi) k iX — -id 

’ 3 


-^(1-4^) 


(Xi,.+XT0-y(Ax)^ 




X’V’-’x.ij - ^X’V.iAx^ = 


2 ^ '^'^3 ’ 


( 4 ) 


( 5 ) 


( 6 ) 


where hij is the tensor-type perturbation generated from the scalar-type perturbation; see equation 21 in jHwang & Noh| ( |2005| l and equation 
126 in jHwang & Noh| ( |2007) ; these works include full third order version of these equations including the vector- and tensor-type pertur¬ 
bations in the context of multiple zero-pressure fluids. In equation l|^, we kept only the second-order tensor-type perturbation generated 
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from the scalar perturbation; this is the only surviving tensor perturbation term in our situation ignoring the linear tensor mode. In order to 
determine Xi, we need equation (|A 6 |l to the second order as 


3A(^) c 

+ C + TC^Xi = ^ 

/ ,2 


(2(^Ax - ^ (‘P.iiX’-’ + 


(7) 


which can be further decomposed as following with scalar and vector contributions X and Yi\ 

H —X i H— Yi, 

a a 


^ 2 A 

^2 


2<pAx- <p,iX’" + {<P,ijX’^ +X,: 


a 


[<P,ijX’^ + X,iA<pJ 
^ + A(px,i - ViA"^V^ (v^.jfeX’'' + XjAipj^ = Wi. 


( 8 ) 


In order to close the system of equations above, we need ip to the linear order. From equation i A5 i to the linear order, we have (with K = 0) 

(9) 


4nGSQ + -«: + c^^</J = 0. 
a a‘‘ 


3 GENERAL RELATIVISTIC ENERGY AND MOMENTUM CONSERVATION EQUATIONS 

By identifying the density and velocity perturbations 

a - -lv7 

Q a a 

we can reveal Newtonian correspondence of our equations derived in Einstein’s gravity. 

We find Xi to second order from equation § as 
c . 

- T 


( 10 ) 


'x=^ (cVx + x‘"') = U + A-1 (VA + Y), 

O 

X = 2ipX ■ u — u • V</9 + - A~^V • (u • VV(p + uAip), 

Y = 2 [u • VVi/p + uA(p — VA~^V • (u • VV(p + uAip)] . 


( 11 ) 


The equations above are valid to the second order with divergence-free vector field Y (that is, V • Y = 0). The linear-order equation reads 

c^X = 

Since 93 appears only with second-order quantities in equations 0,0 and tn), we only need 93 to linear order. From equation 0, we 

find 


2 A . ^ r- Ct 1 

c = —AnGgS -V • u. 

a a 


( 12 ) 


For K — 0 (but with general A), we have yj constant in time; see equation 1 Cl 1 1 for the general solution. In particular for Einstein-de Sitter 
models (K = 0 = A), we have 


<f(= Vv) = -iPx, 


(13) 


where (p„ and are the ifi (diagonal part of Sgip variable in the comoving gauge (v = 0 ) and in the zero-shear gauge (x = 0 ), respectively. 
We have (for general K and A but without anisotropic stress) 


- ^ TT 

‘fix ~ ~ o 


(14) 


where U is perturbed Newtonian gravitational potential (a^ is —()(?oo part in the zero-shear gauge). The metric perturbation variable ip (with¬ 
out taking the gauge condition) completely determines the intrinsic curvature perturbation of our fully non-linearly perturbed cosmological 
space-time, see equation 63 in Noh (2014). 

While the vector contribution is accounted by the presence of Y term via Xi”^ in equation |^, the tensor contribution is accounted by 
hij in equation l|^. The tensor contribution is derived in the Appendix [ b] We find that the tensor contribution to the non-linear evolution of 
density contrast comes through Zij which is given as 

Zij = Nij - 2A-V(iAj'),;, -f J A-" [XXJj -f Si, A) 


Nij = —^ \ u,ijX ■ u 


V ■ u + u ■ — ^Sij [(V • u)^ -f u • (Au)] I 


where A(ij) = \{Aij Aji). 


(15) 
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With these Newtonian identifications, equations 13 and 13 now take the form similar to, respectively, the usual Newtonian energy and 
momentum conservation equations as. 


5 + -V • u + iv • (5u) = - {VS) ■ [2(pu - (VA + Y)] , 


-V ■ 


a 

u H—u 

a 


+ AtvGqS + • (u • Vu) = ^-1 — • V (V ■ u) + 4V 


u ■ Vu-uV • u 

3 


+ -AV u 


u • (VA + Y) - A [u • A"^ (VA + Y)] 


+ 2--u"A-^ . 

a a ^ ' 


( 16 ) 


(17) 


These are the general relativistic energy conservation and (divergence of) momentum conservation equations, respectively, of a zero-pressure 
fluid valid to third order in perturbation. We have assumed a flat background {K = 0) and have taken the comoving gauge; to the linear 


order these equations are known to be also valid on the curved background (that is, in the presence of non-zero spatial curvature A; (Bardeen 

[7^ ). 

Left-hand-sides are the same as Newtonian perturbation equations; the zero-pressure Newtonian perturbation equations have quadratic 
order non-linearity only, thus these parts are valid to fully nonlinear order in the Newtonian context jPeebles||l980^ . Terms in the right- 
hand sides are third order and purely general relativistic corrections in the zero-pressure fluid. Thus, we have exact relativistic/Newtonian 
correspondence to the second order perturbations ( |Noh & Hwang|2004l ). We note that our equations are valid to the third order perturbations 
but in fully relativistic context (e.g., including the super-horizon scale). Compared with our previous presentation in |Hwang & Noh| ( |2006[ l, 
now we have additional Y and Zij terms on the right-hand sides. These correction terms take into account the quadratic-order vector- and 
tensor-type perturbations simultaneously excited by the linear-order scalar-type perturbations. 

We note that while the pure relativistic corrections from the scalar- and vector-type perturbations appear via the linear spatial curvature 
perturbation (p(x, t), the correction from tensor-type perturbation involves linear velocity perturbation u(x, t) only, see equation 1 15 i. 


4 POWER SPECTRA 


We will present the leading non-linear power spectra of the density and velocity field. We introduce the velocity gradient variable ( |Peebles| 

[7^ 

6 = -V ■ u = —K. (18) 

a 

The solutions of 5 and 6 in the Fourier space are presented in the Appendix [C| The density and velocity power spectra are defined as 

{5(ki,f)^(k2,f)) = (27r)®5(®>(ki +k2)P5(fci,f), {6l(ki,f)6((k2,f)) = (27r)®5(^’(ki -f k2)Pe(fci,f). (19) 

We may decompose the non-linear power spectra as P = P 11 -I-P 22 +P 13 , where the subindices ab in Pab{k) stand for two perturbation orders 
out of which the contribution Pab{k) are calculated. Due to the relativistic/Newtonian correspondence up to second order, the relativistic 
Pii (fe) and P 22 {k) terms are the same as the Newtonian ones, and purely relativistic corrections occur from P 13 that involves the third-order 
solutions. We decompose P 13 to the relativistic/Newtonian and pure Einstein parts as P 13 = P^ + P^®, and further decompose P^ to the 
scalar-, vector- and tensor-type contributions as P]^ = Pf^ -f P^^ -f Pm^. 

For density power spectrum, we find 


Pii(fc,f) = \Sdk,t)f 


[ dx 

poo 

/ dr 

7a: -b 3r — lOra:^ 

1-1 

/o 

14(1 -b r^ — 2rx) 


|(5i(fcr, f)|^|(5i(fc-\/l + r'^ — 2rx,t)f, 


Pi 3 {k,t) = ^\5i{k,t)f j 2r(6 - 79r^ -f 50/ - 21r®) -f 3(r^ - l)®(7r^ -I- 2) In 


IL 2 poo 

PMk,t) = -^\Si{k,t)f^ I dr\5dkr,t)f 


' 336r3 


PT{k,t) = £^Mk,t)f’^ 


2r(156 -b 631r^ -b 5l/) -b 3(52 -b 193r^ - 262r'‘ -b 17r®) In ( - - 

\ r + 1 

DO 

Ar\5i{kr, t)\^ 

2r(—36 — 65r^ -b 43r'*) — (36 -b 53r^ — 46r"^ — 43r®) In ( — -^ 

V r + 1 


112r3 


pBV 

03 


pBT 

03 


^ 6(r-b2r®)-b(3-b5r^-b6r'‘)ln^i^q^ 

= ^\3i{k,t)f^^^ / 2r(3-2r^-b3/)-b3(r^-l)^(r^ + l)ln^i^J-J 


( 20 ) 
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Figure 1. The matter density power spectrum with complete relativistic calculation including all scalar, vector, tensor modes up to third order in perturbation. 
We assume no linear-order vector and tensor perturbations, but include the vector and tensor perturbations excited due to non-linear mode coupling and their 
back-reaction on to the scalar modes. We show various components in the matter power spectrum with different colours, and dashed lines indicate that the 
components are negative. The total power spectrum P^{k) is shown as a solid blue line, which lies right on top of the Newtonian non-linear matter power 
spectrum at the same order (P^(fc), green line). It is because the total relativistic corrections black dashed line) are much smaller than the Newtonian 

terms {P22 orange line, P^ (k), pink-dashed line). The relativistic corrections on all scales are dominated by the scalar contribution (P^^ (k) cyan line) 
over the vector {P^^ (k) magenta line) and tensor {P^^{k), yellow line) contributions. 


For velocity power spectrum we find 


Pii{k,t) = \ei{k,t)\\ 


k^ 




pi poo 

/ dx 

/ dr 

/-I 

lo 


r — 7x 6rx^ 

J J Q ^ 14(1 + r'^ — 2rx) 

Pi 3 {k,t) = Ar\&i{kr,t)\^ 


|(5i(fcr, t)\^\Si{k\/l + — 2rx, t)\^ 


168r3 

Pf3{k,t) = - 


2r(6 — 41r^ + 2r^ — 3r®) + 3(r^ — l)^(r^ + 2) In ^ ^ | 

^ Ar\5pkr,t)f 

2r(312 + 527r^ + 207^) + 3(104 + 141r^ - 314/ + 69r®) In ( 'CPll 

\ r + 1 


336r3 


|r- 1| 


PW{k,t) = - — \5pk,t)\ 


> knit)' 


poo 

/ dr\Si{k 
Jo 




' 112r3 


2r(72 + 25/ - 23/) + (72 + / - 50/ - 23/) In 


k-1| 

r + 1 


Pi® 3 ^(fc,t) = J Ar\Si{kr,t)\^ ^^^^^^ ^ 2r(6 + 5r^) + (6 + 3r^ + 5/) In 

Pfi{k,t) = ^\5i{k,t)\^^^^^^ dr\Si{kr,t)\^^^ 2r(3 - 2r^ + 3/) + 3(/- 1)^(1 +/) In 


( 21 ) 


In Figs [^and]^ we present the density and velocity power spectra, respectively, with each contribution represented separately. For the 
cosmological parameters, we use the maximum likelihood parameters in ‘WMAP5+BAO+SN’ column of table 1 of [Komatsu et al.|p009| ), 
and calculate the matter linear power spectrum from the code in|Lewis et al.|([2000ll. Both density and velocity power spectra show that, even 
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Figure 2. Same as Fig.[^ but for the third-order power spectrum of velocity gradient field 9 = (l/a)V ■ u. Again, the relativistic corrections are small for 
all scales, and is dominated by the scalar contribution (k)). 


with the new contributions from induced vector perturbations and tensor perturbations, the pure general relativistic corrections are much 
smaller than the Newtonian non-linear contributions on all scales. 


5 DISCUSSION 

In this work, we presented the third-order perturbation equations (Section 3) with solutions in Fourier space (Appendix C) and the leading 
order non-linear density and velocity power spectra (Section 4) of a zero-pressure fluid in the Friedmann background world model. The main 
new contributions in tbis work are the following. We include the vector- and tensor-type perturbations accompanied by the linear-order scalar- 
type perturbation. As these are simultaneously excited by the scalar-type perturbation, it is necessary to include these for a completeness. 
All these are pure relativistic effects. The full relativistic energy and momentum conservation equations valid to third order perturbation are 
presented in equations and l |17| l. We also present the non-linear power spectra in FigsandThe leading nonlinear contributions to 
density and velocity power spectra via vector- and tensor-type perturbations are smaller than but comparable to the one from pure relativistic 
scalar-type one. We show these pure relativistic corrections are negligible in all scales. Therefore, as a conclusion we proved that overall 
effect from the pure Einstein’s gravity is negligible in the evolution of zero-pressure fluid during the matter dominated era. That is, even 
though the system of equations in general relativity is highly non-linear, there is no significant dynamical contributions arising from the 
non-linearities in general relativity. 

As galaxy surveys in the near future such as DES|^ PF^will be able to measure primordial non-Gaussianities of order /nl — 1, 
it is important to understand the general relativistic effects (that is expected to be important on large scales; therefore, degenerate with the 
primordial non-Gaussianity signal) to better accuracy. While we have shown that dynamically induced non-linear corrections in the matter 
power spectrum are small, there are other general relativistic effects that might be detectable from the large-scale structure observation. 
Eor example, |Bartolo et al.|p005^ studied effects of primordially generated (i.e., inflationary) non-Gaussianity on the later evolution of 
the gravitational potential, and |Bruni et al.| ( f2014| > likewise studied the general relativistic effects through the initial conditions encoded 
during the inflationary era, and similarly considered effects on the non-Gaussianity. Compared with these works emphasizing the second 
order general relativistic effects imprinted from the early universe, our work concerns pure general relativistic third order non-linearity 
generated (dynamically induced) from all three (scalar, vector, tensor) types of perturbations during the matter dominated era. In our case, all 
cosmological parameters are fixed by the A cold dark matter concordance cosmology, and in that case we have shown that the pure general 
relativistic effects on the power spectra turn out to be suppressed in all scales. Another general relativistic effect on the observed galaxy 

^ http://desi.lbl.gov 
^ http://sumire.ipmu.jp 
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power spectrum is the projection effect through the deflection of light coming from galaxies. This effect has been studies extensively to both 
linear order ( |Yoo et al.| 2009[[Yoo| 2010[|Bonvin & Durrer| 2011[|Challinor & Lewis) 20lT||Jeong et al.| 2012{|Yoo| 2014[|Jeong & Schmidt] 

I 2015^ and to second order jBertacca et al.|2014a|b[|Yoo & Zaldarriaga|2014^ . 

There are literatures concerning induced second-order tensor perturbation from the quadratic combinations of linear scalar perturbation 
(Mollerach et al. 2004; Baumann et al. 2007; Ananda et al. 2007; Sarkar et al. 2008; Arroja et al. 2009; Assadullahi & Wands 2009; As- 
sadullahi & Wands 2010; Jedamzik et al. 2010). In this work, we have considered the effects of this induced second order tensor perturbation 
on the third-order scalar perturbation which is demanded to get the leading order non-linear power spectra. We have included the effect of 
similarly induced vector perturbation as well. We note an important difference that we considered the comoving gauge in our calculation, 
whereas all the previous works on the induced tensor perturbation took the zero-shear (often known as longitudinal, Newtonian or Poisson) 
gauge. 

In this work we only considered the scalar-type perturbation in the linear-order, and assumes no linear order vector- and tensor-type 
perturbations. The primordially generated gravitational waves (tensor-type perturbation from inflation, for example) may affect the density 
and velocity power spectra to the non-linear order differently than what has been presented in this work. Although the pure Einstein’s 
gravity effect starts appearing from the third order, the presence of the linear tensor-type perturbation can affect the energy and momentum 
conservation equation from the second order; see equations 39 and 40 in |Hwang & Noh| ( |2005[ l, and equations 125-128 in [Hwang & Noh| 
( [2007^ in the presence of multiple component of zero-pressure fluids; although these equations are valid even in the presence of vector-type 
(rotational) perturbations, as the rotational perturbation only has the decaying mode (in expanding phase) to the linear order, we do not have 
plausible mechanism to sustain the effect of the rotational perturbations. The effect of primordial tensor-type perturbations on the nonlinear 
density and velocity power spectra is an open question left for future investigations. 
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APPENDIX A: FULLY NON-LINEAR PERTURBATION EQUATIONS: REVIEW 

For completeness we present the complete set of fully non-linear scalar- and vector-type perturbation equations in the presence of the 
additional linear tensor-type perturbation. Our metric convention is (Bardeen 1988, Hwang & Noh 2013) 

ds^ = —a^ (1 -I- 2a) (dx^)^ — 2axidx°dx^ + c? [(1 -I- 2<^) Xij + 2hij] dx^dx^, (Al) 

where spatial indices of Xi and hij are raised and lowered by yij as the metric. The metric Xij is the comoving part of the three-space 
metric of the Robertson-Walker spacetime; 6K is the spatial scalar curvature of the background metric Xij (Noh 2014). Here we assume a 
to be a function of time only, and a, ip and Xi are functions of space and time with arbitrary amplitude. We include the transverse-tracefree 
(L^l^- = 0 = hj) tensor-type perturbation only to the linear order; a vertical bar indicates a covariant derivative based on Xij as the metric. 
The spatial part of the metric is simple because we already have taken the spatial gauge condition (to fully non-linear order) without losing 
any generality (Bardeen 1988, Hwang & Noh 2013). The energy-momentum tensor of a fluid in the energy-frame is (Ellis 1971, Ellis 1973) 

Tab — pUaUb 4“ P )Qab 4“ UaUb) , (A2) 


where Ua is the normalized fluid four-vector with u°'Ua = — 1 ; R and p are the covariant energy density and pressure, respectively, with 
tildes indicating covariant quantities. We ignore the anisotropic stress in this work. We introduce (Hwang & Noh 2013) 



7 = 


1 



V^Vf^ 

c^{l+2ip) 


(A3) 


where Vi is an arbitrary function of space and time with the spatial index raised and lowered by Xij as the metric; 7 is the Lorentz factor. The 
complete set of fully non-linear and exact cosmological perturbation equations without taking the temporal gauge condition is the following; 
for derivation, see Hwang & Noh (2013) and Noh (2014). 

Definition of k (perturbed trace of extrinsic curvature): 



a 



1 

A/*(l 4- 2p) 


3^+4 


f 1 

|fc 4- 


x’"T,k \ 

l + 2ipj 


ADM energy constraint: 


3 

2 


3c2 ^ a^{l + 2ifi) 3 ) a a^(l + 2ip)'^ 



A-kG 


(r + R) (7^ 


, \ 3 c pp 

> 2a?{l + 2p)^ 




(A4) 


(A5) 
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ADM momentum constraint: 

2, c ri 

3'^’* aW(l + 2 ip) [2 


(Axi + X^lik) - ^X^lki) + (m + S) ay 


aW(l + 2 :p) 1 V AA 1 + 2 (p 


A/'j ‘Pj 


5 (x"|i + Xi'") - 


<p’-' 


XiP.j + xXi‘P,i + 


(1 + 2(^)2 V'"*’"’" ' 3 

Trace of ADM propagation: 


1 + 2^9 


-3 


1 /a\ 47rG Ac^ 

Af U +:? + !? 




I xV,i + Xip'^ - ^six'°p,k 


1 . „a C^AA/" 


+ —^ K. + 2 — K + 


Af a a'^Af{l + 2 p) 


X + c ^ ^ ^ _ ) +c KjKi. 


a 2 A/'(l + 2 v 9 ) V l + 2 (p 


Tracefree ADM propagation (with linear tensor perturbation): 
;• • 2 A — 2A' fid a 

tlij SHhij — C -;;- hij + ( -77 — + 3 — — K -\- 


cx 


2 yAJ" dt a a'^Af{l + 2 ip) ) \ a^Afll+ 2 ip) 

1 


1 / t I K\ 1 ri fc I “ 

2 (X b- + Xj ) - ^ 5 jX \k - [ X P 

^2 


1 + 2(p 


I \i ^ ci k 

,j 4-Xjy=> - ^OjX ^,k 


a2(l + 2^?) 
SnG 


^ (VV2 4.]a) ^ ^ (VV2 4.]A ) Af 


in + p) 


Lc2(1 + 2v^) 3^^^ f 


aW2(l + 2:p)2 


/ ilk k\i\ 

(^X Xj\k ~ Xfc|iX j 


i I A: k I i 


+ ll2ip (A^^'^XfeVSj - X^'^XjPy + Xmx'^X - XilfeXV''') + _p^2v?) 2 (A‘Xj¥2'"'‘P.fe “ X'XfeV?'V 


b) 


a2(l + 2 ^ 9)2 [1 + 2^3 
Covariant energy conservation: 


V2'V,9 - ^ (¥>'W,2 + V2,2A^'‘ - 


Afu" + -x") V, 

a / 


a 1 

\_dt a(l + 299) 


{^v%k m’‘p,k 1 

a(l + 2 ip) a(l + 2 </ 9)2 ^ 

Covariant momentum conservation: 


/i+(/i + p)'^A/’(3- — K 


^ 1 
dt a(l + 2 ip) 


K + IX) V. 


7 ( = 0. 


1 

ay 


We have 


d 1 

Tvr + 


a(l + 2 p) 

1 [Afc^^ ^ 
fj,-Xp I ay^ 


+ -x") Vfe 

a / 


(ayvi) y-v'^Vi ^ 


a2(l + 2ip) 




72/ a(l + 2 (p) 




Af = ,/1 + 2a + X . = 1 + 5Af, K]Ki = + Xp.) - ^X pX^'p 


1 


1 + 2(p 


a2(l + 2v?) 
i^xV'-’ (xiii +Xiii) - A'pxV., 


aW2(l + 2v9)2 I 2 


1 PIP / \ f i 


+ 


3.V 1*^ ^.9 I T 2y,)2 yX'XiP P,i + gX'xVpVS.P 


(A6) 


(A7) 


(A8) 


(A9) 


(AlO) 


(All) 


As the dimensions we consider 


[gab] = [Ma] = [a] = [yij] = [a] = [p] = [xi] = [xf^l = [wi/c] = Vc] = [ 7 ] = 1, [*“] = [cdt] = [dp] = L, 

\X]=T, [k]=T-\ [v/c]=L, [f^t,] = [Jl] = [gc^] = [p\, [G3=T-^ [A] = [K]=L-\ (A12) 

The perturbation variable tt is a perturbed part of the trace of extrinsic curvature (Kl = —3a/a + k). ADM indicates the Amowitt-Deser- 
Misner (3 + 1) formulation of Einstein’s gravity(Arnowitt, Deser & Misner 1962). Covariant conservation equations are part of the covariant 
(1 + 3) formulation equations (Ehlers 1993, Hawking 1966, Ellis 1971, Ellis 1973). 
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APPENDIX B: TENSOR CONTRIBUTION 

Here we derive the tensor contribution via hij-teim in equation We consider K = 0. We need the tensor-type perturbation generated to 
the second order with the quadratic combination of linear scalar-type perturbation as the source. The tensor perturbation equation generated 
by the scalar and vector perturbation to the fully non-linear order is presented in equation 95 of Hwang & Noh (2013); also see our equation 
l |A 8 ^ ; to the second order, see equation 210 in Noh & Hwang (2004). To the second order in the comoving gauge considering only the scalar 
perturbation, from equation (|A 8 |l, we have 


hij 4 “ 3 " hi-j C ^ hij — C Sij , 
a a? 


(Bl) 


where 


Sij = Hij - 2A -b ^ (V*Vj -b <5ijA) 

1 / ^2 \ 
a^riij = KX,ij - 2X’^X,ijk - 2(pv5,ij - - -5ij ( kAx-“ ‘itpAtp - (p’'^(p,k ) 

(X o V d / 


(B2) 


As we consider hij sourced by quadratic combinations of linear scalar-type perturbation, hij can be regarded as the pure second order. 

In Einstein-de Sitter model, we have a oc oc » 7 ^; 77 is the conformal time with cdt = adij. To the linear order, from equations i A4 1 
l|A 6 ^ and (|A10^, equations l|A 8 |> and l|A10|l, and equation l|A 6 ^, respectively, we have 


n ^ t V 2 ^ 

= 0, - (ax) =‘fi, K = -c -^x- 

n n‘^ 


(B3) 


We can show that the (^-terms contribute as constants (in time) to a^riij (thus to a^Sij as well), and the other terms contribute as a(oc rj^) to 


a?nij\ the growing solution of x is proportional to f(oc a^^^). We change equation iBl ' as 


d'^hij Adhij , , _ 2.-1 


(B4) 


where x = kr} with the wavenumber k defined as A = —k^. We set the source part as 


—a^A ^Sij = x^Aij + Bij, 


(B5) 


where Aij and Bij are constants in time. The homogeneous solution is (Lifshitz 1946) 


, h ■ / \ V / \ sma; cos a; cos a; sinx 

hij oc -ji(a:), -yi{x) oc —5 -^ H- 

a a x'^ x^ x^ x^ 


(B 6 ) 


Using this, the inhomogeneous solution becomes 


hij = [x^ — 10 ) Aij -b Bij = —a^ A ^Sij — IQ Aij. 


(B7) 


Thus, we have 


^X’"^hij =A ^ (a^Sij) =-2H^x’''^^ ^ [a?Sij - k^Bij) , 


(B8) 


where Sij is in equation jB2|; Bij is terms with (p in that equation; we set a?Zij = a^Sij — k^Bij. From these we can show equation Illy. 
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APPENDIX C: MODE-COUPLING ANALYSIS 
Cl Equations in Eourier space 

In Fourier space with /(k) = J equations 1 16 1 and |l7| give 

5(k,t)+ 6'(k,t) = J (fqi j d^g2(5''^^(k-qi2)‘^^i^|^(5(qi,t)6'(q2,t) 

'*’(2^ / / ^^93<5'^Hk-qi23)^^^^^|^5(qi,t)6'(q2,t)v5(q3,t) 

~(2^ / qi2)^5(qi,t)^ [qi • q2X(q2,t)+qi • iY(q2,t)] , 

e(k, t) + 2Jie(k, t) + 47rG^?5(k, t) = f d^q, f d^q25^^'> (k - qi2)^-5^ii-^0(qi. t)e(q2, t) 

j j qfq^ 


(Cl) 


'*’(2^ / 


^(--qr.q2+4 


ai?2 

qi23 • qiqi • q 2 4 


7? 


- gqi23 • q2 ) 6l(qi,t)6'(q2,t)<^(q3,t) 


”'’(2^ / qi2)^6»(qi,t) |?X(q2,t) + ^ [qi • q2X(q2, t) + qi • iY(q2, t)] 

f d?qi f d^q25^^^(k - qi2)%%2JI6»(qi,t)a^Zij(q2,t), 

J J 7i72 


(27r)3 

where equations (U) and give 


(C2) 


Y(k,t) = ^ J J d^q2<5^"'(k-qi2) (^2 

/ ®*"525‘®nk - qi 2 ) 


qi • q 2 3 qi 2 • q 2 qi • q 2 3 qi 2 • qi 52 
5? +2^ gP 2 g2^ qf 


iY(k,t) = 


qi • q 2 , gi qi 2 • q 2 qi • q 2 qi 2 • qi qi 
Q2 2 Ql 2 Q12 2 2 Q12 2 2 

9l 9 i 2 9l <112 <11 


a6»(qi,t)vj(q2,l), 

2 ' 


x2a6'(qi,t)vj(q2,t), 


kkk^ 


1 f kikj 


k,J 


Z^.iKt) = Y,,(k,l) - 2^N^^{k,t) + - +5. 

"«(“•<> /■^«-o-) (jA. - ( 1 + 


q2iq2j 


■Yw(k,t), 


(27r)3 

It is convenient to have 


— 6l(qi,t)6l(q2,t). 


(C3) 


kX(k 


, t) + iY(k, t) = J d^qi J (i^g2(5^®' (k - qi2)a6l(qi, l) 92 (q 2 , t) 

„ /"o qi-q 2 o^(‘ll 2 • q 2 )(qi • q 2 ) + (qi 2 • qi)gi 

‘ 1 - I' - ^ - 2 - 


, o / qi • q 2 , gi 

+ 2g ( q2 —2 -1' qi 

yi yi 


glg^^ij(k,t) =J d^q2 J d^g3(5^®'(k - q23)6'(q2,t)0(q3,t) 


1 + 


q2 ■ q3 
<l2 


i „2 (qi • qs)^ _ ^qi • k /^i 


x^3gr- 


<l3 


fc2 V3 


-qi • k - 


qi-qsk-qsX 1 / (qi • k)^ 2 \ / 1 (k ■ qs) 


ai 




+ gi 


fc^g: 


(C4) 


(C5) 


where we introduced y and z indicating the vector- and tensor-type contributions, respectively; by setting y = 1 = 2 ; we are considering 
these contributions simultaneously excited by the linear scalar-type perturbation. 


C2 Linear solutions 

We set the linear density and velocity fields as 


5i(k,t) = D(l)5i(k), 
ei(k,t) = -FI(l)/(t)73(l)5i(k), 


with / = d In D/d In a satisfying, from equation 1 17 1 
d{HfD) 


dt 


-f 2H‘'fD - A-rGqD = 0, 


or, in a well known form (Lifshitz 1946, Bardeen 1980, Peebles 1980) 
D + 2HD - A-kGqD = 0, 


(C6) 


(C7) 


(C8) 
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where H = a j a. This linear order equation is valid for general K and A in the background (Bardeen 1980), and is also available in the 
Newtonian context (Bonnor 1957, Peebles 1980). Our equation is valid in the comoving gauge and (only to the linear order) the same equation 
is also derived in the synchronous gauge (Lifshitz 1946). Considering general K and A, equation jC8| can be written as (Hwang 1994) 


1 

aa 



(C9) 


with general solution (Heath 1977, section 10 in Peebles 1980) 
S{k,t) = {e - 3K) C{k)H J" J. 


(CIO) 


The transient (in expanding phase) mode proportional to H is absorbed in the lower bound of integration. The coefficient C(k) is normalized 
so that we have (Hwang 1994) 




C(k) 


{l + KH 



(Cll) 


As mentioned below equation l |l2| ip is constant for vanishing K. In Einstein-de Sitter model we have S (X a, a thus the growing mode 
gives D X a. Solutions with general A can be found in Edwards & Heath (1976), Heath (1977), section 13 in Peebles (1980), Weinberg 
(1987), and Heath (1989). 


C 3 Nonlinear solutions 

We expand the higher order density and velocity fields as 


Sn{k,t) = ^2^)3(n-i) / d^9n5^®'(k-qi2...„)74(qi,--- ,q„,f)( 5 i(qi)--- 5 i(q„), 

6 >„(k,f) = J rf^qn5'®^(k-qi2...„)G„(qi,--- ,q„,f) 5 i(qi)---< 5 i(q„), 


(C12) 


with Fn and Gn, the density and velocity kernels, respectively. We have Fi = 1 = Gi. As the pure relativistic contribution appears from 
the third order, the linear and second order solutions are the same as the usual Newtonian ones. 

Now we consider Einstein-de Sitter model. We have D x a, thus f — 1 and 0 = —HS. To the linear order equation \\2) gives 


¥’i(k) 



(47rG£>5 - HO) 


2 


( 5 i(k,f). 


(C13) 
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We then calculate the right hand sides of the energy and the momentum conservation equations, to third order in (5(k) as 

5(k,f)+ 6 '(k,f) = J (fqi j 




+ - 




(? {2n 

5 a'^H^ 1 

2 c2 (27r)' 


^ J d^qi j d^q2 J (k - qi23)(5(qi, f) 0 (q 2 , f)5(q3, f) 

/a./a./ d®g3(5^^'(k - qi23)(5(qi,f)6'(q2,f)(5(q3,f) 


5qi • q 2 


iihs 


<lh<l3 


, o q2 • q3 1.. q23 • q3q2 ■ q3 + q23 • q2?3 


, o ( q2 • qs , gi 

+ 2y ( qi • q3-3-f qi ' q2 ^ 

Q 2 ^2 


g(k, t) + 2H6(k, t) + 47rGg^(k, t) = - J d^qi J (k - qi2)0(qi, f)6'(q2, 

/a./a./ d®g3(5^^'(k - qi23)6'(qi,f)6»(q2,f)(5(q3,f) 


q2qi2 ■ q2 

<ll<l2 


5 1 

'2 c2 (^6 

1 

' 9 9 

mi 


2 , , qi23 • qiqi • q 2 4 

-trqi • q2 + 4 - -qi23 • q2 


91 


5 1 

'2 c2 (27r)6 




/a./a./ d®93(5^^'(k - qi23)d(qi,f)d(q2,f)(5(q3,f) 

, 3 q23 • q3q2 ■ q3 + q23 ■ q2gi ^ , J_ A _ 9123 3 

2 9 I 39 I J 9? V Qh J 


q2 ■ q3 

9i 


,0 qa • qs 1,. qas • q3q2 • q3 + q23 • q2?3 

qi •q23 ( 2- -2 -:^(4y-3) 


91 


9l39i 


+ 2i/ ( qi • q3 ^^ +qi • q 2 % 
92 92 


-2 


/A./<i»®/ d^q3d'®'(k - qi23)6’(qi,f)6'(q2,f)6’(q3,f) 


1 

(27r)® 

'1 2 (qi • qs)^ „qi • q 23 /I 


1 


9? 9^3 


1 + 


q2 ■ q3 


9 I 


7 9i 


9| 


- 2 - 


923 


jqi • q23 


qi • q 3 q 23 • qs 3 , i / (qi • qss) 


9| 


+ 


9^3 


I ^ 

+ 9i 


1 _ (q23 ■ qa)^ 

3 92393 


(C14) 


(C15) 


Now, we use the kernel ansatz in equation l |C12| l to translate the equation above to the equations for the second and third order kernels. 
Using kn = aHjc, we have 


1^+2F2-G2 = ^il^=2l2, 


H dt 
1 dG2 
'HUtr 


91 

_ G2 + ^ (i^2 - G2) = = B2, 

2 9195 




(C16) 




2 --F'2(qi,q2) + 2 

95 953 


J / \ 1 2 

-G2(q2,q3) - k„- 


fn, q2 • q 3 If, 04q 23 • qsqs ■ q 3 + q 23 • q 2 q 3 

"■ ■ “»(,' “ AT “ 5*"’ “ -SI- 

= gl3 + fclfGa, 

1 dGa ^ 47rGq , , 9i23qi ' qzs^ / ^ 

^ ^ =-9!9l^ 


929| ^ <1^30-3 

+ 2y Tqi • q3‘^^^^ + qi • q 2 % 
V 95 95 


5,2 1 [2 
+ 2 "=-^! 3 


2 , , qi23 • qiqi • q 2 4 

— ^qi • q2 + 4- 2 -oqi23 • q2 

o Qi 

p, q2 • q3 , 3 q23 • q3q2 • q3 + q23 • qsqi 

^ 2 4” 2 ^2 „2 


95 


92392 


+ 


91 


953 / 


„ ,0 q2 • qs If, 04 q23 • qsqs ■ qs + q23 • q2q3 

q..q„|2-^--(4!,-3)-- 


+ 2i/ ( qi • qs + qi . q 2 ?| 
92 92 


lXi2 


1 + 


q2 ■ qs 

9 I 


9 i 923 

i„2 (qi • qs)^ _ 2511 • q 23 ft 


391 - 


9| 


9^3 


qi • q 23 - 


qi ■ q 3 q 23 • qs 3 , 1 / (qi • qss) 


9| 


9^3 


+ 9i 


1 _ (q23 • qs)^ 

3 92393 


= B 3 + kfjDs, 


(C17) 
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where ^ 2 , B2, • • • are constants in time. Using k% oc a~^ , the second and third order kernels have the solutions 
F2 = i (5^2 - 2 B 2 ), G2 = ^ (3^12 - 4B2), 

F3 = ^^3-^B3 + ^fc|r(5C3-2i33), G 3 = ^^3 - + yfclr ( 8^3 - 47^3) • (CIS) 

We need symmetrization over wavenumber indices of B2(qi, q2), G2(qi, q2), B3(qi, q2, qs, t) and G3(qi, q2, qs, f). The pure relativistic 
corrections are contained in the coefficients of k% terms. 

This paper has been typeset from a TbX/ KTgX file prepared by the author. 



